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We estimate the perturbative nuclear shadowing in the nuclear structure function F 2 (x,Q 2 ), 
mainly at the kinematic region of HERA using the Glauber-Mueller approach. The contributions of 
the quark and gluon sectors to the nuclear shadowing are estimated. We predict that the nuclear 
shadowing corrections are important and that saturation of the ratio Ri = F 2 (x, Q 2 )/AF%(x, Q 2 ) 
occurs once the shadowing in the gluon sector is considered. 
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^L | The future ultrarelativistic heavy ion collider experiments at the BNL Relativistic Heavy Ion Collider (RHIC) 

qv , and the CERN Large Hadron Collider (LHC) are expected to exhibit new phenomena associated with an ultradense 

,—{ ' environment that may be created in the central collision region of these reactions 0] . The main conclusion emerging 

_ , from the analysis of nucleus-nucleus collisions for RHIC energies and beyond, is that the most of the entropy and 

transverse energy is presumably produced already during very early times (within the first 2 fm after the nuclear 

contact) by frequent, mostly inelastic, semihard gluonic collisions involving typical momentum transfers of only a few 

GeV J2J. However, one is still far from a complete and detailed picture, as reflected by the considerable uncertainty 

in perturbative QCD predictions for global observables in nucleus-nucleus (AA) collisions at collider energies, such as 

particle multiplicities and transverse energy production |3| . The inability to extrapolate accurately from the pp (pp) 

£> ' data to heavy ion AA collisions is due the current poor knowledge about the nuclear and dense medium effects. One 

^sO of these effects is the nuclear shadowing. It is of interest in high energy nucleus-nucleus collisions because it could 

influence significantly the initial conditions in such reactions. Some processes strongly affected by nuclear shadowing 

at collider energies include the minijet production, heavy quarks, their bound states, and dilepton production [0. 

Consequently, the nuclear shadowing is one of the major theoretical issues in modeling the QCD processes in nuclear 

collisions. 

qv . In recent years several experiments have been dedicated to high precision measurements of deep inelastic lepton 

scattering (DIS) off nuclei. Experiments at CERN and Fermilab focus especially on the region of small values of 

the Bjorken variable x = Q 2 jIMv, where Q 2 = — q 2 is the squared four-momentum transfer, v the energy transfer 

and M the nucleon mass. The data [f|||, taken over a wide kinematic range 1CP 5 < x < 0.1 and 0.05 GeV 2 < 

Qh! Q 2 < 100 GeV 2 , show a systematic reduction of the nuclear structure function F^(x, Q 2 )/A with respect to the 

free nucleon structure function F^ (x,Q 2 ). This phenomena is known as the shadowing effect. The analysis of the 

shadowing corrections for the nuclear case in deep inelastic scattering (DIS) has been extensively discussed Q. It is 

motivated by the perspective that in a near future an experimental investigation of the nuclear shadowing at small x 

and Q 2 » 1 GeV 2 using eA scattering could occur at DESY Hadron Electron Ring Accelerator (HERA). 

The deep inelastic scattering off a nucleus is usually interpreted in a frame where the nucleus is going very fast. In 
this case the nuclear shadowing is a result of an overlap in the longitudinal direction of the parton clouds originated 
from different bound nucleons H. Recently, a perturbative approach has been developed to calculate the gluon 
distribution in a nucleus [[7| using perturbative QCD at small x. This approach, known as Glauber-Mueller approach 
is formulated in the target rest frame, takes into account the fluctuations of the hard probe. It includes the shadowing 
corrections (SC) due to parton rescatterings inside the nucleus, and provides the SC to the nuclear gluon distribution 
using the solution of the DGLAP evolution equations [g[ to the nucleon case. As a result the behavior of related 
observables (F/,F^, ...) at high energies can be calculated. 

Our goal in this Letter is to estimate the perturbative nuclear shadowing to the nuclear structure function FA(x, Q 2 ), 
using the Glauber-Mueller approach proposed in M at HERA-A kinematic region. We compute the nuclear shadowing 
corrections to the structure function F% associated with the rescatterings of the qq pair on the nucleons inside a nucleus 
using the Glauber-Mueller approach. Our results are dependent of the behavior of the nucleon gluon distribution. We 
then consider two possibilities. First we contemplate the unshadowed gluon distribution, which we denote the quark 
sector contribution for the nuclear shadowing. Secondly the shadowed gluon distribution, which we denote as the gluon 
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sector contribution for the nuclear shadowing. This contribution is calculated using the Glauber-Mueller approach. 
Our work is motivated by our recent results in the analysis of scaling violations of the proton structure function. The 
ZEUS data for the F% slope presents a turn over which can be successfully described considering the contribution 
of both sectors. Here we estimate the contribution of both sectors to the ratio R\ = F A (x, Q^/AF^x, Q 2 ) and we 
show that the inclusion of the contribution of the gluon sector modify strongly the behavior of this ratio. 

Before we calculate the perturbative SC at F A some comments are in order. As x ~ Q 2 /s, where s is the squared 
CM energy, the data in the region of small x values are for small values of Q 2 (< lGeV 2 ). In this region the 
application of the perturbative QCD cannot be justified. The SC in this region are dominated by soft contributions, 
as can be concluded considering some approaches Jll| [Q (JL3| Q that describe these data. For instance, the 
NMC experimental results Q can be described using the DGLAP evolution equations || with adjusted initial parton 
distributions, following refs. fSBS. Since the initial parton distributions are associated with nonpcrturbative QCD 
this approach intrinsically includes the nonperturbative contributions to the shadowing corrections. Therefore, we 
consider that the nuclear shadowing corrections observed in the NMC 0] and E665 data 0] at small x are dominated 
by soft contributions. The determination of the soft contributions to the nuclear shadowing is not the goal of this 
Letter. 

Let us start from the Glauber-Mueller approach proposed in fij . In the small x region the gluon distribution governs 
the behavior of the observables. In the nucleus rest frame we can consider the interaction between a virtual colorless 
hard probe and the nucleus via a gluon pair (gg) component of the virtual probe. In the region where x << l/2mR 
(R is the size of the target), the gg pair crosses the target with fixed transverse distance r t between the gluons. 
Moreover, at high energies the lifetime of the gg pair may substantially exceeds the nuclear radius. The cross section 
for this process is written as M 



a 9'A = I dz I *Jlw\C?, rt , X , Z )\*o°<*\zrt) , (1) 







where g* is the virtual colorless hard probe with virtuality Q 2 , z is the fraction of energy carried by the gluon and 
ty 9 is the wave function of the transverse polarized gluon in the virtual probe. Furthermore, a 99+ (z, r 2 ) is the cross 
section of the interaction of the gg pair with the nucleus. 

To estimate the nuclear shadowing we have to take into account the rescatterings of the gluon pair inside the 
nucleus. The contributions of the rescatterings were estimated using the Glauber-Mueller approach in ref. |7j. One 
of the main results is that the a 9 A cross section, which can be obtained considering the s-channel unitarity and the 
eikonal model, is given by 

o 9 * A = f dz fil* f^-M* (Q 2 ,r t ,x,z)\ 2 2[l-e-i<^'^ s ^] , (2) 



where bt is the impact parameter and a 9 ^ is the cross section of the interaction of the gg pair with the nucleons inside 
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of the nucleus. Using the relation a 9 (x, Q ) = (An a s /Q )xGa{x, Q ) and the expression of the wavefunction "L 6 
calculated in |16 71], the Glauber-Mueller formula for the nuclear gluon distribution is obtained as 



4 f 1 dx' 
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xG A (x,Q 2 ) = - r=L ^i — 2[ i_ e -i<W,T^]. ( 3 ) 

7T J x X' J 4 TTJ"? J Q -K 

The equation (|3j) allows to estimate the value of the nuclear shadowing corrections in the nuclear gluon distribution. 
In 1 17] the authors demonstrate that a 9 ^ = (3a s (-%)/4) 7r 2 r 2 xGn(x, -4-), where xGn(x,Q 2 ) is the nucleon gluon 

distribution. Therefore, the behavior of the cross section (g) and of the nuclear gluon distribution (||) in the small- a; 
region are driven by the behavior of the nucleon gluon distribution in this region. The Glauber-Mueller formula (0) 
gives us the possibility to calculate the SC using the solution of the DGLAP evolution equation. In this Letter, we 
use the GRV parameterization p9| as a solution of the DGLAP evolution equation. It describes all available nucleon 
experimental data quite well for Q 2 > 1 GeV 2 [|2(J. It should also be stressed here that we disregard how much of SC 
have been taken into account in this parameterization in the form of the initial distributions. The possible SC in the 
initial conditions are nonperturbative contributions to shadowing. 

The use of the Gaussian parameterization for the profile function S(b t ) — (A/irR 2 )exp(—b 2 / R 2 ) simplifies our 
calculations. Then, doing the integral over bt, the master equation is obtained as 

xG A (x,Q 2 ) = ^ / ^- Q ° ±^{C + ln(KG(x',r 2 t ))+E l (K G (x',r 2 t ))} (4) 



where C is the Euler constant, E\ is the exponential function, the function kq(x, rf ) = (3a s A/2R\) it r 2 xGn(x, -4-), 

A is the number of nucleons in a nucleus and R 2 A is the mean nuclear radius. If equation (^) is expanded for small 
kg, the first term (Born term) will correspond to the usual DGLAP equation in the small x region, while the other 
terms will take into account the shadowing corrections. 

The master formula (Q) is obtained in the double logarithmic approximation (DLA) M. As shown in J7|] the DLA 
does not work quite well in the accessible kinematic region (Q 2 > 1 GeV 2 and x > 10~^. A more realistic approach 
must be considered to calculate the nuclear gluon distribution xGa, which implies the subtraction of the Born term 
of (0) and the sum of the GRV parameterization M . This procedure gives 

xG A (x,Q 2 )=xGT ster H. {A)]+AxG G N RV (x,Q 2 )-A^± f f Q ^% x'G% RV {x\Q' 2 ) . (5) 

77 Jx JQl x H 

The above equation implies AxG G J RV (x, Q 2 ) as the initial condition for the gluon distribution and gives 
AxG G ! RV (x,Q 2 ) as the first term of the expansion with respect to kq. Therefore, this equation is an attempt to 
include the full expression for the anomalous dimension for the scattering off each nucleon, while the use of the DLA 
takes into account all SC. In M this procedure was applied to obtain the shadowing corrections to the nuclear gluon 
distribution, demonstrating that the suppression due to the shadowing corrections increases with In 1/x and is much 
bigger than the nucleon case. For Calcium (A = 40) the suppression varies from 4% for In 1/x = 3 to 25% for 
lnl/x = 10. 

The expression (|^) estimates the SC to the gluon distribution using the Glauber-Mueller approach. The modification 
in the nuclear gluon distribution represents the contribution of the gluon sector for the nuclear shadowing. 

The Glauber-Mueller approach can be extended to the nuclear structure function considering that in the rest frame 
of the target the virtual photon decays into a quark-antiquark (qq) pair long before the interaction with the target. 
The qq pair subsequently interacts with the target. In the small x region, where x « l/2mR, the qq pair crosses 
the target with fixed transverse distance rt between the quark and the antiquark. Considering the s-channel unitarity 
and the cikonal model and following the same steps used in the case of the nuclear gluon distribution, the nuclear 
structure function can be written as 

9 n f 1 t9 

F 2 A (x,Q 2 ) = ^fJ24 f^ d -I±{C + ln{K q {x,ri))+E x {K q {xy t ))}, (6) 

3tH ^r> J i 7rrf 
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where n q — 4/9 kq. Similarly as made in the case of the nuclear gluon distribution, to obtain a more realistic approach 
the Born term should be subtracted and the GRV parameterization should be added. Therefore the nuclear structure 
function is given by 

F 2 A (x, Q 2 ) = F 2 A (x, Q 2 )[eq. (6)] - F A (x, Q 2 )[Born] + AF 2 (x, Q 2 )[GRV] , (7) 

where F A (x , Q 2 )[Bom] is the first term in the expansion in n q of the equation (ft), and F2(x,Q 2 )[GRV] — 
J2 U d s e q [ x l( x ' Q 2 ) + xq(x, Q 2 )] + F 2 (x, Q 2 ) is calculated using the GRV parameterization. The charm component of 
the structure function is calculated considering the charm production via boson-gluon fusion 119]. In this Letter we 
assume m c = 1.5 GeV. 

The expression (R) estimates the SC to the nuclear structure function using the Glauber-Mueller approach. We see 
that the behavior of F 2 is directly associated with the behavior of the gluon distribution used as input, which is usually 
described by a parameterization of the parton distributions Jl9| . In this case the shadowing in the gluon distribution 
is not included explicitly. We denote the quark sector contribution for the nuclear shadowing the effect in the nuclear 
structure function obtained using the Glauber-Mueller approach and a usual parameterization (unshadowed) as input. 

The total modification in the structure function associated with the nuclear shadowing is the sum of the shadowing 
corrections in both sectors. The main difference between quark and gluon sectors stems from the much larger cross 
section a 9 ^ — (9/4)cr|^, i. e. kg — (9/4)k 9 , which in turn leads to a much larger gluon shadowing. Therefore, at small 
values of x the contribution of the gluon sector cannot be disregarded in a reliable calculus of the nuclear structure 
function. In a general case, the shadowing corrections for F 2 should be estimated considering also these corrections 
for the gluon distribution, i.e. in the quark and gluon sectors. Recently, we propose a procedure to estimate the SC 
in both sectors to describe the turn over in the F 2 slope observed in the ZEUS data J2lJ . Using the procedure here, 
we estimate the nuclear structure function considering both sectors (quark+gluon sectors) of NS using the solution 
of (H) for the nucleon case (^4 = 1) as an input in the expression (M). 

Using the expression (|7j) we can estimate the nuclear shadowing corrections. We assume that the mean radius of the 
nucleus R\ is equal to R\ = 2/5 R%/g, where Rws is the size of the nucleus in the Woods-Saxon parameterization. 
We choose Rws — r^As , with r = 1.3 fm. We calculate the ratio 



Rl AFi(x,Q*) ' (8) 

where F%(x, Q 2 ) is the structure function for a nucleus with A nucleons and F^ix, Q 2 ) is the proton structure function. 

In figure ffl we present our results for the ratio R\ as a function of In 1/x at different virtualities and considering the 
different sectors for the Calcium case (A = 40). In Fig. [l](a) we present the results of the NS for the ratio considering 
only the quark sector, while in Fig. [y(b) we present the results of NS corrections considering the quark+gluon sectors. 
We can see that in both cases the perturbative nuclear shadowing is an important correction to the behavior of 
the nuclear structure function. The suppression varies from 1% for In 1/x = 3 to 20% for In 1/x = 10. Therefore 
our result demonstrates that the suppression due to shadowing corrections is proportionally smaller in the nuclear 
structure function than in the nuclear gluon distribution. This result is in disagreement with the result obtained in 
ref. p2[ which uses the Gribov formula and a Regge parameterization for the diffractive structure function. We believe 
that this disagreement occurs since that approach is dominated by soft physics. Furthermore, we can see that when 
only the contribution of the quark sector is considered no saturation in the ratio R\ is observed, which agrees with 
the result obtained in |2J|. However, if the contribution of the gluon sector is considered, i.e. if a shadowed gluon 
distribution is used, the ratio saturates in the HERA-A kinematical region. This is a new result, since that in general 
the saturation in ratio R\ has been characterized as essentially a nonperturbative effect [g3[. Therefore, a saturation 
in the perturbative region (Q 2 > 1 GeV 2 ) evidentiates the large shadowing corrections in the gluon sector. 

In figure @ we present our results for the ratio Ri as a function of Q 2 at different values of x and different sectors 
for the Calcium case (A = 40). As the behavior of the ratio i?i in the region Q 2 < 2.5 GeV 2 is strongly dependent 
of the initial conditions p3], we present only our predictions in the region of larger Q 2 . We can see that the ratio is 
dependent on Q 2 in the region Q 2 < 10 GeV 2 and that at large Q 2 the nuclear shadowing corrections diminish but 
do not vanish, i.e. it is not a higher twist effect. In addiction, for x > 10~ 3 the curves associated with the quark and 
quark+gluon sector coincides, which represents that gluon sector do not contribute in this kinematical region. The 
dependence of the ratio i?i on Q 2 occurs mainly at small values of x, but the general behavior of the curve do not 
depend on the sector considered. 

In figure we present our results of the ratio i?i as a function of A at different values of Q 2 and x for the different 
sectors. We can see that the nuclear shadowing corrections increases with the growth of A and with the decrease of 
x. The nuclear shadowing corrections for the ratio R\ are larger for smaller Q 2 values and smaller if the gluon sector 
is considered. 

In summary, the perturbative nuclear shadowing corrections for HERA-A are analyzed in this Letter in the context of 
the Glauber-Mueller approach. This approach predicts important results for the perturbative nuclear shadowing: the 
nuclear shadowing corrections are very large in the HERA-A kinematical region; saturation of the nuclear shadowing 
corrections will occur in this kinematical region if the contribution of the gluon shadowing is considered; the shadowing 
in the nuclear gluon distribution is larger than in the nuclear structure function. We expect that these results motivates 
the running of nucleus at HERA in the future, since the precise determination of nuclear shadowing corrections is 
fundamental to estimate the cross sections of the processes with nucleus which will be studied in the future accelerators 
RHIC and LHC. 
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FIG. 1. The ratio Ri = F 2 (x, Q 2 ) /A F% (x, Q 2 ) as a function of Inl/x at different virtualities and different sectors: (a) 
quark sector and (b) quark+gluon sector. In the Calcium case A — 40. 
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FIG. 2. The ratio Ri = F^(x, Q 2 )/AF%(x, Q 2 ) as a function of Q 2 at different values of x and different sectors (quark and 
quark+gluon). In the Calcium case A = 40. 
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FIG. 3. The ratio Ri = F 2 (x,Q 2 )/AF2(x,Q 2 ) as a function of A at different values of Q 2 , different sectors (quark and 
quark+gluon) and (a) x — 1CP 3 , (b) x — 1CP 4 . 



